Regularization of multi-soliton form factors in sine-Gordon model 



O 

(N 

\0' 



p. 

a 

(N 
> 

00 

o 



T. Palmai* 

Department of Theoretical Physics, 
Budapest University of Technology and Economics, 
H-1111 Budafoki ut 8, Hungary 

April 17, 2012 



Abstract 

A general and systematic regularization is developed for the exact solitonic form factors of exponential 
operators in the (l+l)-dimensional sine-Gordon model by analytical continuation of their integral represen- 
tations. The procedure is implemented in Mathematica. Test results are shown for four- and six-soliton 
form factors. 



1 Introduction 



Form factors (matrix elements of local operators) are important quantities in quantum field theories. It is a 
O ^ remarkable feature of certain two-dimensional field theories (integrable models), that their S-matrices can be 
obtained exactly in the framework of factorized scattering theory [I] [2]. Furthermore, in integrable models there 
is a rather restrictive set of equations satisfied by the form factors (that is the form factor axioms I31I2]), which 
Ci • makes it possible in many cases to obtain them exactly as well. For instance, in the sine-Gordon model all 
form factors of exponential operators are known O HI [5]. The spectrum of the sine-Gordon model consists of a 
soliton-antisoliton doublet and their bound states, called "breathers". While the breather-breather form factors 
can be given explicitly (see e.g. [3]), the solitonic ones, in general, are only known in terms of some highly 
^ ■ non-trivial integral representations. In addition, the integrals converge in a limited domain of the parameters. 
\^ , In this paper we give a regularization procedure to calculate the solitonic form factors in the sine-Gordon model 
00 ■ for arbitrary choice of the parameters. The regularized multi-soliton form factors could them be used to obtain 
[ correlation functions of direct physical interest, e.g. in condensed matter physics [S]. 

The outline of the paper is as follows. In Section 2 the sine-Gordon model along with its exact form factors 
are reviewed and based on [3] integral representations for the form factors of exponential operators are given. 
Section 3 is devoted to the analysis of a certain function which appears in the integral representations. Giving 
this function's asymptotic series and identifying its poles make it possible to analytically continue the integral 
representations. In Section 4 explicit formulae are provided for the four-soliton form factors. Section 5 is devoted 
to discussion of test results, while Section 6 is left for conclusions and outlook. 

><: 

cs; 2 Form factors in the sine-Gordon model 

2.1 Definitions, S-matrix and the form factor axioms 

The sine-Gordon model is defined by the classical Lagrangian 

ia^^av + ^cos(/3^) 



L = I dx 

Define the parameter 



(2.1) 



which is relevant in the low-energy description of the theory. The spectrum of the quantum theory contains 
the soliton-antisoliton doublet and their bound states, the "breathers". The number of breather states [Bi, B2, 



., Bjv) is bounded, there are N 



of them. For our purposes it is enough to consider only the solitonic 



particles of the spectrum, indexed in the following with ± (soliton-antisoliton). 
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2.1.1 S-matrix 

The S-matrix for the sohton-antisohton doublet reads 

/ S 



St Sr 
Sr St 



(2.3) 



S J 



with the non-zero elements [T] 

S+Z (6) 
S+t (6) 

where 



= 511 (6) = 5 (9), 
= szl (6) = St (6) = 5 (e) - 



sinh Y 



sinh 



st+ (e) = 5fl(e)-5(e) 



sinh^ 



sinh 



iTT-e ' 



(2.4) 
(2.5) 

(2.6) 



5(e) = -(-i)~n 



fc=l 



X exp 



y sm(et)- 



2sinhlii2^e-^<* 



2 sinh 4^ cosh ^ 



(2.7) 



which is independent of the integer N , however the integral converges in a larger domain of C for iV > 0. 

2.1.2 Form factors 

Consider the matrix elements 

(©i' ■ • ■ out {A'- (e:„) . . . A"^ (e'l) (Gi) . . . (e„))^„ (2.8) 

of the local, hermitian operator O between asymptotic states. The form factors are defined by 

^ai...a„ (61, . . . e„) = {Q\0\Aa, (61) . . . (e„)),„, (2.9) 

as the matrix elements of the operator between the vacuum and an n-particle state. Crossing symmetry implies 

^l\zlz (©i' • • • e;jei, . . . e„) = J-a,...a„(-bi)...(-6™) (©i^ • • ■ ©'i + ©™ + i^) , (2.10) 

which is understood as an analytic continuation. The form factors can be reconstructed based on the following 
axioms. 

1. Analyticity and the physical matrix elements. J^qi...q„ (8i, . . . Qn) is analytic in the variables Qi — Qj 
inside the physical strip < Im O < 27r except for simple poles. It is the physical matrix element when 
all @i are real and ordered as 81 < 82 < ... < 8„. 

2. Relativistic invariance. The form factors satisfy 

^a,...a„ (81 + Z, . . . e„ +Z)= e^^<-"^Ta,...a„ (81, . . . 8„) , (2.11) 

where S (O) is the spin of the operator O. 

3. Watson's theorem. The following symmetry properties are satisfied 

J^ai...ajaj + i...a„ (81, . . . 8j, 8j + i, . . . 8„) = 5^^^^^+'. (8j+i - &j)Tai...CjCj + i...a„ (81, . . . 8j + i, 8j, . . . 8„) , 

(2.12) 

.Fa,...a„ (81, . . . 8„ + 2TTi) = e^""^^^*) J",,.,, , (8„, 81, . . . 8„_i) , (2.13) 

where the latter is understood as an analytic continuation and uj (O, ^) is the mutual non-locality index 
of the operator O and the "elementary" field, upon which the whole operator product algebra can be 
constructed. 
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4. Kinematical poles. J^ai...a„ (©i, • • • ©«) has simple poles at 0„ = Oj + itt with residues 

^iJ'a[...a'....a'^ (^©li ■ • ■ i • ■ • ©n-l^ |<5aJ • • ■ <5aj_ J S'o"_ Jci (©n-1 — ©j) 'S'a"_2C2 (©n-2 ~ ©j") • ■ ■ 

x4:^ar-^(©,+i-©,) 
x5:-;!rM©.-©.-i)4:j---fc}- 

(2.14) 

In the absence of bound states these are the only singularities of J-ai...a„ (©i, ■••©«) in the strip < 
Im©j < 27r for real 

2.2 Integral representations of multi-soliton form factors 

In [3] it is proposed that the 2n-particle form factors of the exponential operator e^"""^ in the sine-Gordon model 
can be represented by 



^ai...a2„(©l, • ■ • , ©2n) = (0|e"^'^|^.,„ (©2„) ■■■A,, (©i)) 

= Ga{{Z^,„ (02„)...Z,, (©l)» =^a-F.i... „(©!,..., 02n) (2.15) 



where X]i=i = because of charge conservation and Qa is the vacuum expectation value of the exponential 
operator [H [7j . The operators Z± (©) are defined by 



Z+(0) = ,/^^e^e^^(^), 



(2.16) 



Z_(©)^,/.^e-^ e^^y l^eV---. 



^ii-^-iT)(7-e) ,0(0) 



-e 



^g(l-f-iT)(7-e) ,0(9)^^0(7) 



27r 



(2.17) 



Since <?!)(©) and are free fields the averaging ((. . .)) is performed by the multiplicative Wick's theorem, 
using 



((e"^(e2)gi0(ei)^^ = G(©i -©2), 
((e'0(e2)g-#(ei)^^ = 1^(61 - ©2) = 

((e-»^(02)e-'^(®i))) = G(0i - ©2) = 

The appearing functions and constants are as follows. 

G(©) =iCi sinh (^j exp ■ 



w{e) = - 



cosh (©) ^'^'^ I jQ 



1 



G(0i-©2-f)G(0i-©2 + f)' 

1 

W (01 - ©2 - f ) (©1 - ©2 + f ) ■ 

dt sinh^ t (1 - ^) sinht(C - 1) 1 
1 sinh(2t) cosh(t) sinh(^t) J ' 

dt sinh^ t (1 - ^) sinh t (C - 1) 1 



sinh 2t sinh 



G(©) = -|?sinh(^^±^^) sinh(©), 



(2.18) 
(2.19) 

(2.20) 



(2.21) 

(2.22) 
(2.23) 



^ dt sinh^(|)sinht(e-l) I 

^ /o * sinh(2t) cosh(i) sinh(^t) J 



G(-*7r), 



C2 = exp 4 



sinh^ (I) sinht(C- 1) 
1 sinh(2t) sinh(^t) 



Ty(f)esin(|) 



(2.24) 
(2.25) 
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The integration contours appearing in the consequent expressions for the form factors are such that the 
"principal poles" of the M^-functions are always between the contour and the real line. (We define the "principal 
pole" of W (9) as the pole located at 9 = -^). 

For the two-particle form factors it is only necessary to evaluate Eq. (|2.15p for two Z operators. Let 

A = - (^^ + , then the result is 

{{Z+ iQ2)Z- (9i))) = !^et(e2-eOG(ei - 92)6-^^^ 

X / ^e'^'^^(7-02)M^(9i-7)-e"i5^ J ^e^'^M^ (7 - 93) (7 - 9i)| , (2.26) 

(92)^+ (9i))) = ^ef (ei-e^)G(9i - 92)e-^^^ 

X |e^5^ j ^e^^W^ (02 --f)W (9i - 7) - j '^e^'^W (7 - 92) W (9i - 7)| , (2.27) 

The four-particle form factors can also be obtained through evaluating Eq. (|2.15p with the result 



((Z,, (94)^^3 (03)^.. (92) ^.,(9i))) = ^^^^G(93-94)G(92 -94)G(9i-94) 

XG (92 - 93) G (9i - 93) G (9i - 92) Ja,<y,<y,a,. (2.28) 

where 

7 - piT.'l=l'^^'S'^„~AY.a,=-l^iT ('2 29^ 

and /(Ti(T2 0-3fT4's are given by 





= e^P(/22,/3l) 


— P (/22, ^4o) 


-P(/3i,/3l)^ 


h e" 


"^^P(/3i,/40), 


(2.30) 




= e^P(/i3,/3i) 


-P(/22,/3l) 


-P(/l3,/40)^ 


h 


"^P(/22,/40), 


(2.31) 




= e^P(/04,/3l) 


-P(/l3,/3l) 


-P(/04,/40)^ 


h e" 


-|^P(/l3,/40), 


(2.32) 


/++__ 


= e^P{Io4,Il3) 


-P(/l3,/i3) 


-P(/04,/22)^ 


h 


"I^P(/l3,/22), 


(2.33) 




= e^P(/04,/22) 


-P(/l3,/22) 


^ P {IqA:, ^ 


h 


-|^P(/l3,/3l), 


(2.34) 




= e^P(/l3,/22) 


— P (/22, ^22) 


-P(/l3,/3l)^ 


h 


'^P{h2,hl). 


(2.35) 



The integrals 7^^ have four components, lij^k fc = 1, . . . , 4 and the operation P is defined by 

P (a, 6) = e~ (ai&i — 0262) — e~~ (0363 — 0464) . (2.36) 



/ij^fe's read 



(2.37) 



/o4,fc = y e(^+"'=)^Ty (94 - a;) (93 - x)W [Q2 - x)W (9i - a;) dx, 

/i3,fc = j e(^+"'=)^iy (x - 94) W (93 - x) (92 - x) (9i - x) dx, (2.38) 

/22,fc = f pJ-^+^'-^^Ty (x - 94) Wix- 93) (92 -x)W (9i - x) dx, (2.39) 

g(A+a,)x^ (a; - 94) W{x^ 93) (a; - 93) W (9i - a;) dx, (2.40) 



i31.fc 



/40_fe = / e('^+"'=)^W^ (a; - 94) T4^ (a; - 93) W {x ~ &2) W {x ~ 9i) dx (2.41) 



with «! = — 1 — i, a2 = 1 — |-, 013 = — 1 + i, (14 = +l + i coming from writing G{x) as the sum of four 
exponentials; the contours are as before. 
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In general, the 2n-particle form factor is realized as 



71 



X / < 



n (7- e,+„) J] W (sign (j - z) (7, - 6,)) I . [] G (7, - 7,) 



(2.42) 



"("-1) 



and the sum of 4 2 exponentials containing 7i's. 



The last product gives the numerical factor 

All in all, we have (n + 1) combinations of the M^-functions, which must be integrated over with some exponential 
factors. Note that the exponential factors do not alter the structure (e.g. the poles) of the integrands. The 
other kinds of 2n-particle form factors can be obtained e.g. through the symmetry properties of form factors 
(Watson's theorem). 

The problem with such integrals is that they diverge for either 



or 



Rea > — 
2 



1 ^ 
Rea <-- + -. 



(2.43) 



(2.44) 



For such choices of a the integrands have essential singularities at Re x ^ ±00. In the next section we prove that 
W (x) has an asymptotic series in exponentials of x, therefore the divergent integrals can always be analytically 
continued to obtain a finite result. Our strategy is to first deform the integration contours to the real line, then 
extract the divergent terms of the integrands in the form of exponentials and give their contributions exactly 
by the analytic continuation rules 



1 



a e 



+ 00 



a e 



exp (ax) dx = 

exp (ax) dx = 

Then if the integral is expected to be analytic and f (x) dx exists we have 

f[x)dx^Y.-+ (/ (^) - [/] (^)) + / /(^)^^ 

00 j J -00 Jo 

for some / (t) admitting an asymptotic expansion in exponentials, 

/(x)= ^ a,e"'" + 0(e"+")= [/](x)+0(e"+"), x ^ -00, a+ > 0. 

a,<0 



(2.45) 
(2.46) 

(2.47) 
(2.48) 



The previous equation defines the function [/](x). The case when f (x) dx = 00 is similar. 

It should be noted that for some combination of the parameters, the analytic continuation may still produce 
an infinite result, that is in the case a^ = for some i. This happens e.g. for the integral /22.1 when ^ = \- 
These infinities, however, must and indeed do cancel out from our end results, the form factors, therefore the 
a^ = terms in the asymptotic series should be omitted before making the analytic continuation prescribed in 

3 Analysis of the VK-function 
3.1 Asymptotic series 

The function W (x) is given by (|2.22p . The asymptotic series of cosh(x)~-'^ reads as 



cosh(x) ^ = 2e 



l-e 



Re X 



(3.1) 
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where s=signRe x was introduced for convenience. The exponent of the remaining part of W{x) can be rewritten 



as 

[°° dt sinh(C-l)t / , , f2tx\ . , N • f'i.txW 

Differentiate the previous formula with respect to x and obtain 

2 r"" , sinh(^-l)i / , , . (2tx\ . , , /2te\\ , 
c?^ . , . \ . , . cosh (2t) sm + « sinh (2i) cos = (3.3) 



TT Jo sinh (^i) sinh (2t) \ \ tt J \ n 

The asymptotic series of the Fourier integrals was first obtained by |8j, where an analog to Watson's lemma for 
Laplace transforms was discussed. Given a function q (t) with the asymptotic series near t — 0: 

oo 

g(t) = ^6„r+^-i (3.5) 

n=0 

with some < A < 1, the asymptotic series of 

/•CO 

F{x)= q (t) e±**^rfi (3.6) 
Jo 

as a; ^ oo is given by 

oo 

= ^6„e±^("+^)r(n + A)a;-"-^ + 0(e~''^) , x ^ oo, (3.7) 

where O (e~'^^) denotes corrections "beyond all orders", i.e. exponentially small terms (/i > 0). Applying this 
construction to the integrals occurring in Eq. p.3p we get 

2 ^^ sinh (^ - 1) t ^^^^ ^.^ C^)^ -i^Im lime^^r(A) + 0(e-''i-) = i^ + 0(e-''i-) (3.8) 
TT io sinh(^t) ^ ^ V 7^ / 2C A^o ^ ^ ^ 2^ ^ ^ ^ ^ 



2 r,,Binh(^ l)t^^ /2^A ^ ^3_g^ 



and 

TT Jq sinh(^i) \ TT J 

considering that q (t) is odd in the first integral and even in the second which implies that all but the first term 
of the first integral disappears because of the Im/Re operation, respectively. For simplicity it was assumed, that 
X is real. It is easy to check our statements remain true if this condition is relaxed. 

To find the exponentially small terms the integral is evaluated by the residue theorem which yields only the 
residues times 27ri since the integrands are of small enough order on the half circles = {z : ±Im z > 0, |2:| = R} 
as i? -H- cx) (where is associated with Re x > while with Re x < 0) . The result is 

i_l + g cot + ,-(2fe+l).. ^ g ^-2ks. 

^ k=l ^ ^ k=l 

2^ /'7Tk\ f2TTk\ 2. 

cot P^~'' 



.-gsm^^Jcot^^)e-T^^ 



|^(-l)^+^sin((i^il^) e-f^^ 

(3.10) 

When integrated with respect to x one gets the exponent of the VF-function: 



e-1 
2^ 



-sx ■ 



fc=i 



TTk 

T 



cot 



27rfc 



IS 



e « 



E 

fc=i 



1 fi'^+icot , for odd k 



'ksx 



k 1 7^ 



for even k 



The integration constant, Cs is determined from the relation 



Gix) 



1 



W{x+f)W{x-f 
and the explicit form of the function G(x), that is 



(3.11) 



6 



(5(a;) = -^Csinh ) sinh(e) , 



with the result 



- 4e^^ = 



4i 



2^ 



e 25 . 



Now we are ready to give the asymptotic expansion of W (x): 



(3.12) 



(3.13) 



W{x) 



4i 



e 2« e 25 



EM 



-2lsx 



n 

k=l 



00 / 



^ Zl 

./=0 



^ Zl 



where the coefficients depend only on ^ and s=sign (Re x) and are expressed as 



and 




IS sm 



Trk 



for odd A:, 
for even k. 



Rex s ■ 00, 

(3.14) 

(3.15) 
(3.16) 



Note, that, strictly speaking, our expansion is Umited to the case of irrational ^ parameters since otherwise 
the coefficients ak and bk always become infinite for some k. However, such infinities can be shown to cancel 
out. 

Let ^ = where ni and n2 are relative primes. We have singular afe's whenever = odd, which 
immediately implies ni = even and n2 = odd and the A''th singular a-coefficient is indexed by where N 
is necessarily odd. On the other hand bi is singular if ^ = even while I is odd, implying again rii = even 
and 712 = odd and the A^th singular 6-coefficient is indexed by Nn2- Now in both cases the A^th singular term 
contribute terms of order e"^"^^ to the exponent of W (x). All that remains is to show that the A''th diverging 
coefficients cancel each other. To see this let ^ = ^ (1 + s) or equivalently 



resulting in 



m ^ ni(l + e), (3.17) 
na n2(l-e), (3.18) 



1 cos(7r7Vn2) 1 1 

Keajvn,/2 = cos(7riVn2/2) ~^ ~ Nm cos (7riVn2 (1 - e) /2) ' ^ ' 



■'Nn2 



JV(ni+n2) + l 

1 ,jv„,+iC0s(^jVni/2) 1 (-1) 



Nn2 sin(7riVni/2) iVn2 sin (TriVm (1 + e) /2) 

which if expanded in e yield 



(3.20) 



Reajvni/2 = &JVn. +0(£) = (-1)™ -r4 - + 0{e), (3.21) 

which agrees for £ — >■ 0. Because of the cancellation the following rules can be formulated for rational ^'s: 

is Ik 

afc = - — for — = odd, (3.22) 
k ^ 

6fc = for A; = odd and fc^ = even. (3.23) 

3.2 Poles 

With the asymptotic expansion at hand the divergences of the integral representations can be readily remedied. 
However there is another issue with the integrals containing W-functions. In fact, W{x) has a number of poles 
on the line Re a; = 0. When the integration contour is fixed (which is the desired scenario), poles can cross it 
and one needs to analytically continue the result by adding the residue contributions of the crossing poles. In 



7 



the followings we determine the poles of W{x). The poles of W{x) are easily extracted from the identity ([3], 
but also follows from a similar representation of G{x), given in |10) ) 



^ r (1 + H^-|±^ j r M + — 



cosh a; , 2k-^\^ 



k=l 



( /■°°dte-'^^*sinh2i(l-^)sinhi 1 

X exp <^ -2 / . ; ^ 7: , ) ■ 3.24 

1 Jq t sinh2ismh^t J ^ ' 

They originate from the poles of the gamma functions and the roots of coshx. It is apparent that Rex — for 
every pole. We are interested in the poles of 14^ (x — xq) which cross the real line (or the original integration 
contour) when ^ is decreased. The previous equation yields two infinite series of poles given by 

e + 2fc-^ + ^^ - (3.25) 
^ + 2fc - i - ^ (3.26) 

2 TT 

n being a non-negative integer and k being a positive number. We have the following estimates for the series of 
poles: 

Imxi,fe,„ > (!^_^, (3.27) 

Imx2,fe,„ < ^ (^-28) 

For |Imxo| < fvr the only poles that can cross the real line are 

Xn = Xi^l,n=iTT(^i-^ ., 71=1,2,.... (3.29) 

and no poles can cross the original contour, which intersects the Imx — line at x — xq = — — e, e — ^ 0. With 
reference to the form factors, note that because of Watson's theorem it is enough to give a calculation method 
when all the rapidities satisfy |Im0i| < tt. Thus it is not necessary to analyze further the ^-dependent poles of 
W {x — xq), that is covering the case |Imxo| > tt. 



4 Implementation of the four-soliton form factor formula 

As an example of the machinery outlined in Section 2 in this section the implementation of the four-soliton 
form factors is discussed. Implementing Eq. (|2.28p is non-trivial only in the calculation of the integrals Uj.k- 

First, by Cauchy's theorem we deform the integration contour to the real line. For this we need to identify 
the poles between the real line and the original contour, which consists of the the principal poles and (possi- 
bly) several ^-dependent poles (given by Eq. (|3.29p ) of M^-functions. The principal poles give the following 
contributions to l22.k'- 

p,,k = -^'^^^^^"^^^^^"^^)^(0i-04 + f)M/(er-e3 + f)Ty(e.-er-|), (4.i) 
P2,k - -^e(^+-)(«^+^)M^(^e2- 64 + ^)14^(62- 63 + y)M^(ei-e2-f), (4.2) 

P3,k = -^e(^+"'=)(«-^)T4^ (^63 - 64 -y) 1^(^62 -63 + ^)14^ (^61- 63 + y), (4.3) 
P,,k = -^e(^+-)(«-^)M/(64-63-f)vK(62-64 + y)M/(6i-64 + f), (4.4) 
if Im6i > — f , Im62 > — f , Im63 < Im64 < +-|, respectively. The ^-dependent poles yield 
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Ni 


-64 


-xn)w{ei 


-03 


-x„)WiQ2 


-01 


+ x„) 




N2 


-64 


-Xn)Wie2 


-03 


-x„)W{Qi 


-02 


+ Xn) 




71=1 


-64 


+ xn)wie2 


-03 


-xn)wiei 


-O3 






N4 
71=1 


-63 


+ Xn)Wie2 


-04 


-xn)wiei 


-O4 





where 7Vi,2 = L(ImOi,2/7r + 1/2) /CJ and ^¥3,4 = [(-Im 03,4/7r + 1/2) /^J , x„ are defined by Eq. and r„ 

is the residue of W (x) at x„, calculated numerically by the definition: 



rn = lim {x~Xn)W{x) 

X^Xn 



(4.9) 



Second, the integrals, 7^^^ with the deformed contours are evaluated by the analytic continuation formula 
(|2.47p . In conclusion one gets 

/O ^00 4 

{A{x)^[A]{x))dx+ A{x)dx + Y,{Pn.k + X^,k), (4.10) 



with 



A (x) = exp [(A + ak) x]W{x- O4) (x - O3) W (O2 - x) (0i - x) 
= [yl](a;) + 0(e"+^) = ^ a,e"'^ + 0(e"+^), Rex -cx). 



(4.11) 
(4.12) 



Qi<0 



Upon generalization to the 2n-soliton form factors the only non-trivial component of this procedure is the 
determination of the residues picked up when deforming the contour. For an integrand 



we have 



with the residue contributions 



^(x) = e^"[]W(s,(x-0,)), s,;=±l 

i=l 

/oo 
A{x)dx + P + X, 
-00 



ImO; + - 



0i - 0j - 



(4.13) 



(4.14) 



(4.15) 



and 



N Ni 



n — 2silm 0i 



27r^ 



(4.16) 



X = ^ ^ 27rzr„e^(«'+^'''") J] W {s^ [0, - 0, + s,r„]) , iV, 

i— 1 71=1 j^i 

To conclude this section, we give some details of the Mathematica [TT] package SGFF . M. After the above, only 
one element remains that is not straightforward in the implementation: the calculation of the asymptotic series 
of the product of several M^-functions. Simple products of the series quickly produce an intractable number of 
terms, most of which are inaccurate (higher order terms, to which further orders in the series of the constituent 
functions would contribute). Our solution makes use of 2-by-n matrices containing the coefficients and the 
exponents of the terms in the asymptotic series. We have the following key procedures in the package SGFF.M. 

• PRO [2t, calculates the 2-by-m matrix corresponding to the product of asymptotic series 2t, *8 including 
orders only with accurate coefficients. 

• SHI [21, a] generates the 2-by-n matrix corresponding to the asymptotic series of / (x + a) from that of 
fix) (i.e. 21). 
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• AInt [21] gives the integral of the asymptotic series corresponding to 2t on the negative half-hne. 

• Asyf un [21, x] yields the value of the asymptotic series corresponding to 21 at x. 
In terms of these procedures /22,fc is calculated as 

rO 



'22. k 



AInt [21] 



(A (x) - Asyf un [21, x]) dx + / A{x)dx + ^ (P„.fe + X„^k) , (4.17) 
-°° ■'0 n=l 



with A{x) as before, and 

21 = PRO [PRO (PRO [SHI [W* , -64] , SHI [2IJ*, -63]] , SHI (<2U, -83]] , SHI [2IT, -61]] , (4.18) 

2If being the asymptotic series of W{x) for Rex > 0. The remaining integrals are performed by the routine 
NIntegrate. 

From a practical point of view one should note, that the second term in (j4.17p can be numerically unstable. 
On the other hand, provided that the truncated asymptotic series is a good enough approximation of A{x) for 
a; < the contribution of this term can be neglected altogether. Therefore, we omitted this term from our code 
and supposed that the input rapidities are big enough for this to cause no harm. This can be assumed safely, 
since Lorentz invariance implies that the rapidities can be shifted by an arbitrary real number. 

Considering now the general case of 2n-particle form factors, we give the asymptotic series of 

N 

A{x)^e^-l[W{s,{x-e,)), (4.19) 

diverging for a: — >■ —00, as 



W , = +1 
2n*, s, = -1 



21 = PRO.^i [SHI [Co,. (25) , -9,]] , Co,. (2U) = 

The main functions available in SGFF.M are to calculate the two-, four- and six-soliton form factors. 

• FF2[0i, 82] gives the two-particle form factors 

{F+_(ei,, - e2,.),i^-+(ei,, - 62,.)} 

where 81 and 82 are arrays of the same length with elements lda,i (a = 1,2, « = 1,2,...A^) rapidities 
where the two-soliton form factors are to be evaluated. 

• FF4[8i, 82, 83, 84] gives the four-particle form factors 

{F__++(8i,,;, 82,., 83,,, 84,0, ^-+-+(. • •), ^^+— +(• • ■),F-++-{- ■ F+-+-{- • •), ^++— (• • •)}• 



• FF6[8i, 82, 83, 84, 85, 85] gives the six-particle form factors 



(81^,, 82,., 83, 




©5,17 ©6,4, 


),^^++-+— (■• 


•), ^^++- 


-+-(• 


■■).F++ 


— +(■ 






F+_+_+_(... 




f+-(..-) 


, F+_+__+(.. 


^^+— 


^+-(.. 


•),^+- 


-+-+(•■ 






F+— ++{... 


), + - 


f+-(---) 


, J^-++-+-(•• 


^^-++- 


-+(•• 


•),^-+ 


-++-(•■ 






F-+-++{... 


), ^^- + - 


-+-+(■••) 


, F__+++_(.. 


.), 


^-+(.. 


.),F__ 


+-++(•■ 




-+++(•••)} 



• FF2p[8i,82] gives the two-particle form factors for physical rapidities, i.e. ones with imaginary parts of 
±7r. 

• FF4p[8i, 82, 83, 84] gives the four-particle form factors for physical rapidities. 

• FF6p[8i, 82, 83, 84, 85, 86] gives the six-particle form factors for physical rapidities. 

Note that accurate results can only be expected when all rapidities have big enough positive real parts and 
imaginary part in the interval [— 7r,7r]. The functions calculating form factors only at physical rapidities are 
considerably faster compared to the general ones if the form factors are needed in more than one points; they 
calculate the necessary W^-function values for the integrals only once as part of the initialization. 

The parameters that can be specified in SGFF.M are the following, which can be edited e.g. in Mathematica 
before loading the package. 
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• ^ is the IR parameter (|2.2p . 

• aover/3 is the ratio of the parameter a appearing in the operator O — e*""^ and the UV parameter /3 of 
the Lagrangian. 

• NN is the regularization parameter for the G~ and W^— functions denoted by N in the formula p.24p . 

• Na is the maximum number of terms treated in the individual asymptotic series in the formula (I3.14p . 

• Ni is the number of interpolation points used to calculate the integrands of type (j4.19p . When evaluating 
the form factors at general rapidities, mainly Ni determines the time of evaluation. However, it is this 
parameter that influences the accuracy the most, as well. A safe choice is Ni — 2000. 

• £ is a technical parameter for the calculation of the residues (14. 9p . e — x — Xn- 

• aa and bb are the lower and upper bounds of the integrals of the type A[x)dx in (|4.17p 



5 Tests 

The four-particle form factor F |_+ was checked against the free fermion point result (omitting the vacuum 

expectation value) 

{{Z+ (84) Z+ (83) Z_ (62) Z_ (81))) - sin^ (V2na>j eV2a(e4+e3-e.-ei) ^^ -^^ 

sinh {^^^) sinh (^^^) 
^ cosh i^'^) cosh i^'^) cosh (^^^) cosh (^^^7' 

We do not show test results for this formula since our calculations agreed with the exact results to the machine 
precision (of 15 digits). 

Also, we investigated whether the numerically obtained form factors satisfy the form factor axioms. In the 
four-particle case the equation 

Fa,a2a,a^ (©1 + 83 + Z, 83 + Z, 84 + 2:) = F^,^,„,„, (81, 82, 83, 84) (5.2) 

must hold, which was checked. Watson's theorem is another axiom, e.g. in the form 

(81, 82, 83, 84) = S+t (83 - 82) (81, 83, 82, 84) + S+Z (83 - 82) (81, 83, 82, 84) , 

(5.3) 

and 

F^^++ (81, 82, 83, 84 + 2m) = e2"-i^+__+ (84, 81, 82, 83) , (5.4) 

uj = ^ being the mutual non-locality index. 

The residues of the kinematic poles of the four-particle form factors were also checked by: 

i ^ lim (84 - 82 - in) F__++ (81, 82, 83, 84) = F_+ (81 - 83) [S+Z (83 - 82) - e^^'^ SZZ (82 - 81)] . 

(5.5) 

Testing the kinematical poles is especially important: for the cases when the two-soliton form factor is known 
explicitly (e.g. for half-integer ^), equation (|5.5p gives the only check that is independent of numerical integrals 
and their analytic continuations. E.g. for ^ = 1 the two-particle form factors are known to be |i4j 

G(-i7r) \2J ^sinhf^J 

In Tables 1 and 2 we listed test results for the four-particle form factors. One can see that magnitude of the 
error varies greatly for different scenarios. This is because we work in fixed precision (double precision) and the 
integrals appearing in the formulas can assume values of very different magnitudes and rounding errors can get 
magnified. 

We also implemented the 6-particle form factors. The numerical results agreed with the exact results for 
the free fermion point, where the 6-particle form factor F reads 

-I sin" {V2na) eV^°(Q«+Q5+e4^e3-e.-ei) ^-^^^ {^^) sinh (^^) 

©4—01 ^ / 64 — 62 \ „„„V. I 64 — 63 ^ / 65—61 \ / 65—62 \ 



cosh ^!4_^ cosh cosh cosh cosh 



2 ) \ 2 ) \ 2 ) \ 2 ) 

sinh (^i^) sinh (^^2^) sinh (^2^) sinh (^i^) 
■ cosh (^^) cosh (^^) cosh (^^) cosh (^^) ^^■'^^ 



11 





LHS 


RHS 


(Iir^. e = 2.23 
(|5.2p. 4 = 0.34 
((Ol). e = 2.23 
(lOl). ? = 0.34 
(ISa. e = 2.23 
CT. e = 0.34 


0.45330- 1.4092i 
0.00089 — 0.05U 
0.453360- 1.4093198i 
0.0009063 - 0.04937438^ 
-0.04255089122137 + 0.03246926430660^ 
-0.043292833089+ 0.00219194033^ 


0.45336- 1.4093i 
0.00091 — 0.049i 
0.453358- 1.4093196i 
0.0009065 - 0.04937441i 
-0.04255089122139 + 0.03246926430663^ 
-0.043292833083+ 0.00219194037i 



Table 1: Comparison of the LHS's and RHS's of the form factor axioms (|5.2p (where z = 1 was taken), (15. 3p . 
(|5.4p in the four-soHton case. The rapidities were chosen to be 9i — 7.6, 82 — 7, 83 = 7.2, and 84 6 — in. 
In ah the tests % — j was set. 





LHS 


RHS 


(15^,^ = 2.23 
(l^.f = 0.34 


0.8211182 + 0.7147548i 
-0.2812726 + 0.0213804i 
-0.4726029- 0.6620907i 


0.8211175 + 0.7147545? 
-0.2812724 +0.0213801? 
-0.4726070 - 0.6620917i 



Table 2: Comparison of residues of four-particle form factors with exact results (|5.5p . We took ^ = 1 and for 
the rapidities 81 = 7.6, 82 = 7, 83 = 7.2, and 84 = 7 + 10"*^ + in. 

In addition, the following kinematic pole equation was tested: 

I lim (86-83 -z^)F___+++(8i,82,83,84,85,86) = 

(81, 82, 84, 85) [S+Z (85 - 83) S+Z (84 - 83) - e^""^" (83 - 81) SZZ (83 - 82)] . 

(5.8) 

Watson's theorem requires e.g. 

F—+++ (81, 82, 83, 84, 85, 86) = S+Z (84 - 83) +-++ (81, 82, 84, 83, 85, 86) 

(84 - 83) +++ (81, 82, 84, 83, 85, 86) . (5.9) 

Test results for the six-particle form factors are listed in Table 3. 

It interesting to note, that the use of multi-soliton form factors extends to the calculation of soliton-breather 
and breather-breather form factors by virtue of the bound state pole axiom [3]. In case of higher breather- 
breather form factors, using soliton-antisoliton form factors can be preferable: e.g. to calculate the Bn-Bm 
form factor one needs to evaluate either the (n + m)-Bi- or the four-soliton form factors. 

6 Conclusions and outlook 

We established a method to obtain the multi-soliton form factors numerically in the (l+l)-dimensional sine- 
Gordon model. The form factors are known in terms of integral representations, whose domains of convergence 
were extended by analytical continuation. In order to do this we needed the asymptotic series of the TV-function. 
Detailed formulae were only shown for the four-soliton form factors, however the number of treatable particles 
is not limited by the procedure. Test results obtained by the code provided for the four- and six-soliton form 
factors were shown. 

Based on the formalism developed for finite volume form factors in [T^ [T3] , a program to investigate the 
sine-Gordon form factors is currently underway [TU], which can now be extended to multi-soliton states [TS]. The 
present formalism is also expected to be relevant to boundary form factors |16j , for which finite size corrections 
have been developed in [IZ] and applied to sine-Gordon theory in [TB] . Future applications will also include the 
calculation of finite temperature correlation functions based on the formalism developed in \W[ [201 ?IT\ . 





LHS 


RHS 


dEll), e = 1-17 
(1511. C = 0.34 
(I5a),e = 1.17 
(OD. e = 0.34 


-0.50782662 - 2.333030973i 
-0.3945330 - 0.3095434i 

-2.84279- 1.63925i 
-0.0151096- 0.147483i 


-0.50782660 - 2.33030977i 
-0.3945333 - 0.309543H 

-2.84263- 1.63902i 
-0.0151107- 0.147442i 



Table 3: Comparison of the LHS's and RHS's of the form factor axioms in the six-soliton case. For (|5.9p the 
rapidities were chosen to be 81 — 2.1, 82 — 1.9, 83 = 6, 84 = 5.9, 85 = 1.2, 86 = 5.5 + in, while for (|5.8p . 
81 = 2.1, 82 = 1.9, 83 = 5.9, 84 = 1.2, 85 = 6, 86 = 5.90001 + in was taken. In all the tests | = | was set. 
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